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Abstract 

We present a pedagogic derivation of the electromagnetic field established in a dielectric material 
by an impinging external field. We consider the problem from the point of view of the physical 
mechanism involved at the microscopic level. The internal field emerges when the material is 
thought of as an assembly of atoms in vacuum, each of them being polarized by the external 
incident field and by the re-radiated fields of all the other polarized atoms of the material. In this 
way, each atom becomes itself a source of secondary radiation that adds and interferes with all the 
other internal fields (including the internal extension of the externally impinging field), contributing 
to build up the total internal field within the dielectric material as well as the externally scattered 
field. 

This picture naturally leads to a connection between the microscopic properties of the material 
and its index of refraction, that describes the dielectric response to the applied electromagnetic 
field. Calculations also show that the incident radiation is extinguished inside the dielectric by an 
equal but opposite field generated by the total induced polarization currents, and is substituted by 
a macroscopic electric field which propagates with the new speed c/n. 
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I. MAXWELL’S EQUATIONS AND ELECTROMAGNETIC WAVES IN DIELECTRICS. 


Maxwell’s equations 1 provide a compact and elegant derivation of why a light beam, 
entering a transparent dielectric material, propagates inside it with a reduced speed. 

Using Maxwell’s equations (in gaussian units), it is straightforward to derive the wave 
equations for the electric field E(r,t) in vacuum (identical equations hold for the magnetic 
field B(f, t )): 


V 2 E- 


1 d 2 E 
c 2 dt 2 


0, 


( 1 ) 


characterised by dispersiveless propagation with constant speed c, the light speed in vacuum. 

When, instead, an electromagnetic wave propagates in a dielectric material (we think of 
it as homogeneous, isotropic, and non magnetic, and in addition that it responds linearly 
to the incident radiation), a polarization current is induced given by Jp = dP/dt, which 
makes Eq.Q inhomogeneous through the appearance of a source term dJp/dt = d 2 P/dt 2 \ 


1 d 2 E _ 4vr d 2 P 
c 2 dt 2 c 2 dt 2 


( 2 ) 


In Eq.([2j), E(f,t) and P(r,t) are macroscopic fields of Maxwell’s equations in dielectric 
materials, connected by the constitutive relation 


P = XE 




( 3 ) 


where x and n are the frequency dependent electrical susceptibility and index of refraction 
of the material, which also are macroscopic parameters of the theory. Using the second of 
the above relations, Eq.([2]) becomes 

o -* 1 d 2 E (n 2 - 1) d 2 E 

c 2 dt 2 c 2 dt 2 

The source term on the right hand side of Eq.Q reveals an interesting form since it splits 
into two terms. A first one is relative to a propagation with the light speed c in vacuum and 
it cancels the corresponding term on the left hand side. A second one re-establishes, after 
that cancellation, the homogeneous structure of the wave equation in the material with the 
field that propagates with a reduced (and dispersive) speed c/n. 
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II. A MICROSCOPIC APPROACH OF LIGHT PROPAGATION IN DIELECTRICS. 


Wave propagation in its fully general manifestation is an extremely complex phenomenon 
for every type of wave motion we consider. An exceptionally lucid and lively presentation 
of this complexity is given by Feynman in an interview by Christopher Sykes that is posted 
on YouTube 2 . This complexity invites a further investigation into wave phenomena, taking 
different points of view in order to understand and clarify them at a deeper level. 

For example, the previous analysis of light propagation in dielectrics does little to clarify 
the physical mechanism involved at the microscopic level. For this deeper explanation, the 
principle that we must consider when a light beam enters a dielectric is the atomic nature of 
matter. A piece of dielectric material is a volume of vacuum space filled with a collection of 
atoms, each polarized and behaving like a charged harmonic oscillator 3 (an oscillating dipole 
moment) when acted on by the external incident radiation and by the fields re-radiated by 
all the other oscillating dipoles of the collection. 

This picture, equivalent to the treatment with Maxwell differential equations as far as 
the final result is concerned, is mathematically very different since the atomic description 
of matter in interaction with radiation naturally leads to a representation through integral 
equations. From their solutions we get the expression of the index of refraction n of the 
material (expressed in terms of the microscopic parameters of the atoms, of their density, 
and of the frequency of the incident field). In addition, the solutions straightforwardly 
include an equal but opposite field, generated by the total induced polarization currents, 
that extinguishes the incident held inside the material (a fact known in optics as extinction 
theorenW^), and its substitution with a new internal held propagating at a reduced speed c/n , 
despite the fact that the incident beam and the re-radiated waves by the atoms propagate 
with the vacuum light speed c in the interatomic vacuum within the material. 

We shall present this point of view, basing our treatment on purely classical electrody¬ 
namics, but closely following the neat and concise presentation given by M.L. Goldberger 
and K.M. Watson in their technical monograph on collision theory®. 

Let us consider a dielectric medium (it may be glass) made of N impenetrable identical 
atoms, kept together by their bonding interactions and distributed randomly within a volume 
V of arbitrary shape with the only restriction that the radius of curvature of its boundary 
at each point must be large compared to the wavelength A of the incident electromagnetic 
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radiation. We further consider N to be a large number and the atomic density p to be 
approximately constant and equal to N/V within V. 

The medium is placed in the path of a plane electromagnetic wave with wave vector k, 
whose wavelength A = 2n/\k\ is much greater than the atomic dimensions (for example, for 
visible light we have that the wavelength is ~ 10 4 the atomic dimensions). We also require 
that the radiation frequency be significantly less than the frequency corresponding to the 
lowest excitation energy of the atoms (cu << uj a p = (E a — Ep)/h), so that the interaction of 
the incident radiation with the atoms may be considered as essentially elastic. The response 
of the dielectric material to the impinging radiation is then associated with the scattering of 
light by the atoms, and we anticipate that the electromagnetic propagation through it will 
be described by its index of refraction. 

The microscopic, local electric held (monochromatic of frequency oj, with complex time 
dependence e ~ lU3t ) is the real part of E n (r n ) and has the general structure of the multiple 
scattering equations 

N 

Enij'n) Eq {t n) T ^ ' E radm {vf ) (5) 

m(^n)=l 

where E 0 (f n ) is the incident (complex) electric held, and E radm (f n ) is the held radiated by 
the m-th atom, all these helds being evaluated at the position f n of the n-th atom in V. 

This is a complex formula, with the self-consistent signature of multiple scattering prob¬ 
lems, since the electric interaction of the m-th atom at the n-th atom position, described by 
the second term on the right hand side of Eq.([5]), is itself dependent on the held acting on 
the m-th atom, which has exactly the same form of Eq. (J5]) . 

Having assumed that the incident radiation wavelength is >> than the atomic dimensions, 
the dipole approximation applies in this case, and in the second term in Eq.(|5]) we may use 
the classical scattering formula for E radtJn f n ), given by 8 
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where p m is the electric dipole moment of the m-th atom induced by the local held E m (f m ), 
f nm = r n — r m , and V n is the gradient operator acting only on the n-th atom coordinates. 

With the assumptions of isotropic dielectric and linear response, the induced electric 
dipole moment of an atom is connected microscopically with the local electric held acting 
on it by the scalar atomic polarizability, 
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A simple classical formula for the atomic polarizability can be found using the classical 
charged harmonic oscillator model 7 . An easy calculation leads to the result (for a single 
resonance frequency) 

e 2 /m 
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(cap — cu 2 ) 2 + 7 2 w 2 

where e, m,cuo ,7 are the oscillator characteristic parameters and oo is the frequency of the 
electric driving force. 

We can now go over to a continuum representation of matter, replacing the sum over m 
in Eq. (|5]) with an integral over the volume V of the dielectric except for the exclusion of a 
sphere So with atomic radius and centered at the n-th atom. Since we made the assumptions 
that atoms are impenetrable, we are guaranteed that only the chosen atom is excluded as 
required by the sum in Eq. (|5]). We thus obtain (in the following, r n and r m should be 
regarded as continuous vectorial variables) 


A'n(^’n) A 0 (r n) T P I d T rn E r ad,m(j’n) 

JV-S o 


(9) 


Eq. (|9]), together with expressions ([6]) and ([7]), represents the integral equation for the 
field E n . Its solution will give the electric field at every point inside V. To do this, we 
proceed according to the arguments of classical macroscopic electrodynamics, defining the 
macroscopic field in the medium, E(r), as 


E(r) = E 0 (r) + p / d 3 r m E radjm (r), 

JV-cyl 


( 10 ) 


where in this case the integral is over the volume V excluding a very small cylindrical, 
needle-shaped region which includes the point Eand which is parallel to E(r). 

Subtracting (10) from ([9]), and using again this last equation to eliminate the local held 
E n (r n ), we obtain 

p d T m E rac i J m(p'n) T/(r n ) E 0 (r n ) p d T rn E ra d^m{j‘n)- (H) 

Jv-So JSo-cyl 

We thus see that the complex expression of the dipolar interaction among the atoms on the 

left hand side, that also appears in Eq.(|9]), is equal to the right hand side combination of 

the incident held, the internal macroscopic held, and a term that reproduces the one on the 
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left hand side except that the integral now extends to the (atomic-sized) volume contained 
between the sphere S 0 and the needle-shaped surface cyl that intersects it. 

To evaluate this integral, we make a further simplification based on the fact that the 
wavelength of the incident field is >> than the atomic dimensions. This allows us to neglect 
the term ik in Eq.Q as compared with l/r nm (a quasi-static approximation), and to set 
e lkrnm = 1. We can further consider p m = aE m (f m ) independent of r m within So, thus 
replacing p m with p n . Then Eq.(|6|) becomes for small r nm : 

Pm X r nm 
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( 12 ) 


In the final expression of Eq. (12), we have retained only the term that gives a non-zero 


contribution to the integral on the right hand side of Eq. © (see Appendix). So the 
volume integral can be carried out, using standard vector analysis results, transforming it 
to a surface integral. This vanishes on the needle-shaped surface of the cylinder, and only 
receives a contribution from the integration on the surface Sq. The result is 

Pn ' f nm ~\ d/T 


d ^’m^rad,mij’n) 


da r 


~Pn 


(13) 


J So-cyl J So-cyl 

Remembering the definition of the polarization vector and using Eq. (|7j) for the electric 
dipole moment, we can write P{f n ) = pp n = paE n (r n ), and this can be used with the 
equations ( 0 , and (p~ 3 | to obtain the final result 


P{f) = paE 0 (r ) + pa 


-> _ 47T -> 

E(r) - E 0 (r) + —P(r) 


(14) 


This is the self-consistent solution of the integral equation (J9]) we were looking for, which 
evidences how the microscopic response of the material translates into a macroscopic polar¬ 
ization, the source term for all that propagates within the dielectric and also emerges out of 
it (internal field, scattered field, as well as, for special geometries with flat interfaces of the 
dielectric material, refracted, reflected, and transmitted fields). It is also worth noting at 
this point that no hypotheses have been made about the directions of the incident oscillating 
electric field and of its wave vector, so that the obtained relation is general. 

When we only consider the internal fields, it is immediately evident from © that one 
of the effects of the atomic dipolar re-radiations (the square bracket in this equation) is to 
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cancel the incident radiation field from within the dielectric: this is the expression of the 
extinction theorem. 


In addition, Eq. (14) shows the content of multiple scattering of the atomic dipolar re- 


radiations as they reproduce the macroscopic polarization also on the right hand side of it, 
in case of dense materials. After the cancellation of the incident field, and considering that 
within an isotropic material the polarization and the electric field vectors are parallel, we 
obtain 

pa „ n 2 


P = 


-E = 


1 — ^fpa 47 t 


-E 


(15) 


where for the second expression the definition of index of refraction has been used. This is 
the promised connection between the microscopic parameters and the index of refraction. 


In case of low density materials (such as gases), we have pa << 1, and Eq. (15) simplifies, 
providing a simpler relation between the atomic polarizability and the index of refraction. 
The equation becomes n 2 — 1 = 47T pa and, since (n — 1) << 1 in the given approximation, 
it assumes the final form 

n — 1 + ‘In pa. (16) 

In this case the enhancement effect due to multiple scattering is absent. 


Since Eq. (16) is a very often used expression for the refractive index in the elementary 


expositions, we would like to add few more comments on it. In the theory of scattering, 
when the above approximation applies, it is common to describe the scattering off a scatterer 
(like an atom or, its classical counterpart, a charged harmonic oscillator) by the scattering 
amplitude f(k,k'), where the incoming and outgoing wave vectors, in the case of elastic 
scattering that we considered, satisfy the relation \k\ = \k'\ = k. Consider a dielectric 
material in the form of a very extended but thin slab, with thickness L « 1/k, placed in 
the plane z — 0. From what already assumed, the atoms within it have a very low packaging 
density (pa << 1 is satisfied). On this slab a plane electromagnetic wave (assume unit 
amplitude, for simplicity) is incident along the direction z. The electric field computed in a 
point z » L assumes the form of a superposition of the incident field and of the scattered 
fields in z by all the atoms of the slab': 


E(z) = e ikz + e ikz2 ^^f(k, k) 
k 


~ e ikz e i^f(k,k) 


( 17 ) 


where in the final results appears the forward (zero angle) scattering amplitude f(k, k) = 
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/( 0 ), and the expression of the final formula derives from the the assumed approximation. 
But the exponent of the last exponential of Eq.(Jl7|) has a very familiar expression in ele¬ 
mentary optics, being equal to i(n — 1 )kL, so that, by comparison, we get the well-known 


expression 


n = 1 + 


(18) 


and, lately, comparing this expression with Eq.(16), we obtain the relation between forward 
scattering amplitude and atomic polarizability: 


/(o) = k 2 a. 


(19) 


There exists a remarkable relation (of general validity for all the scattering processes, not 
only for light), called the optical theorem , that connects the forward scattering amplitude 
and the total extinction (scattering plus absorption) cross section. For the case considered 
of elastic scattering, this result is easily derived. Having chosen the incident radiation to be 
of unit amplitude, the measured intensity at a distant forward point z is proportional to 


|£( z )| 2 = le'^rO* 0 !! 2 = e -^/m/(o) = e~‘” L 


( 20 ) 


where the last expression only includes the scattering cross section, since scattering is the 
only loss mechanism that we considered during the crossing of the slab along the propaga¬ 
tion direction of the beam. We could as well describe intensity reduction in the forward 
direction as interference effect between the incoming beam and the scattered beam. It is 


now immediate, comparing the exponents of the two last expressions in Eq.(20), to get the 
optical theorem in this case: 

o' = y/m/(0). (21) 


and from it, using equations (19) and (|8j), to obtain the scattering cross section from a 
charged oscillating dipole: 

,4 


fitTrr^ 

a = Anka" = 0 
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3 (ufi-wty + 'YW (22) 

In Eq.(22), r 0 = e 2 /me 2 is the classical electron radius, and we used the relation® 7 = 
2e 2 u 2 /3mc 3 for the radiation damping of Eq. (| 8 ]) . This important result, describing the 
Rayleigh scattering, incidentally explains, at an elementary level, the blue color of the sky 
we see when our line of sight forms a right angle with the white sunlight crossing the 
atmosphere, as well as its red color at sunrise and sunset. 








In the opposite case of dense material, the relation between the atomic polarizability 
and the index of refraction is more complex and the formula of Lorentz-Lorenz is obtained 


solving Eq.(15) for the polarizability: 


n 2 — 1 
n 2 + 2 


47T 


(23) 


III. NOTES ABOUT THE EXTINCTION THEOREM AND THE FIELDS GEN¬ 
ERATED BY THE POLARIZATION CURRENTS 

The previous results, in spite of the fact that they are well-known, are only seldom 
discussed in introductory electromagnetism and optics courses. This, notwithstanding the 
fact that some of them could easily be deduced from simple experimental observations. 

In Fig. (|TJ) an elementary experiment of optics is shown, that is commonly done to demon¬ 
strate Snell’s laws of reflection and refraction. A light beam is incident at an angle on a 
plane air-glass interface, and reflection and refraction beams are produced, all in the same 
plane of incidence. Using his exceptional physical intuition, Feynmarf® interprets this sim- 



FIG. 1. Experiment of a light beam striking obliquely an air-glass interface. Incident, reflected, 
and refracted beams are shown. The dotted line would be the direction of the incident wave in 
glass if it could go straight into it. 
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pie experimental situation in the light of the atomic model of matter to deduce, with the 
aid of the energy conservation law, the characteristics of the reflected beam. 

He notices that a beam along the dotted line in fig. (|Tj) , which continues the incident beam 
straight into the dielectric, is not observed. It should, only if the atoms of the glass were 
frozen in their positions and could not oscillate under the driving action of the incoming 
beam, so that the beam can propagate along the dotted line into the interatomic vacuum 
with speed c. But the atoms do polarize and oscillate driven by the incoming primary held 
and by the secondary re-radiated fields by the multiple scattering atomic actions. The total 
induced polarization current must thus generate, in addition to the refracted beam, another 
vacuum beam that exactly cancels the incoming beam inside the dielectric. This is the way 
Feynman uses the extinction theorem (without ever naming it), just running the argument 
backwards: a vacuum beam inside is not observed, so, by superposition, there must be an 
internally generated vacuum beam that extinguishes it. 

Equally interesting is to study the role of the polarization direction (perpendicular or 
parallel to the plane of incidence) of the incident beam in this experiment. As is well- 
known, these two electric held polarizations lead to different expressions for the reflection 
(and refraction) coefficients. Feynman approaches the problem of deducing the different 
contributions for the two cases according to the idea that only the helds radiated by the 
projections of the oscillating dipoles perpendicular to the propagation directions of both 
reflection and extinction helds are non-zero. He carries out the main calculation steps, also 
offering hints to arrive at the hnal answers that he writes down. A detailed derivation 
along these lines of thought, with some further interesting extensions, have recently been 
published 11 . 

Finally, there remains one more question to be noted about the experiment in hg.Q. The 
standard derivation for the reflection and transmission coefficients, for both perpendicular 
and parallel polarizations, is by matching boundary conditions at the air-glass interface for 
the macroscopic helds of Maxwell’s equations, and of course that works very well except 
that it does not clarify an important point. In fact, the reflected and transmitted helds 
appear, in that derivation, to only depend on what happens at the interface, a fact that 
seems to be also supported by examining the experimental situation in hg.Q. However, 
this is not what we understand with our microscopic calculation, which suggests instead 
a different physical interpretation. What generates all of the internal and external helds 


10 


(except the incident one) is the total polarization current that is set up into the dielectric 
material. If the material has a plane interface, as in hg.Q, then the fields radiated by all 
the atomic dipoles in V, coherently organise themselves in order to generate (in addition to 
the extinction held) fields that propagate along the refraction the reflected directions. This 
sounds more likely as a bulk effect rather than a surface effect. 


Appendix A: Derivation of the results of (12) and (13). 


The following standard vector identities involving the operator V hold: 

V x (Ax B) = A(V ■ B) - B(V ■ A) + (B ■ V)A - (A - V)B 

V(A • B) = (A ■ V)B) + (B ■ V)A + A x (V x B) + B x (V x A) 

If A is a constant vector, then the terms that apply derivatives to A are zero, obtaining 

Vx(IxB) = A(V ■ B) — (A- V)B 
X7(A-B) = (A- V)B + Ax (V x B). 

The term (A ■ V)B can be eliminated by the two previous expressions, getting 
V x (A x B) = - V(A -B)+Ax (V x B) + A(V • B). 


The complete expression for the Eq.(12) is thus: 
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The volume integral of Eq. ( fTTj ) is thus composed by the three terms that appear in Eq.(Al), 
and it can be evaluated using the geometry shown in Fig. (J2]) . The integral of the first term 


of (Al) is converted to a surface integral. The integral over the needle-shaped surface of the 


cylinder gives a negligible contribution, while the one on the surface S'o results in the con¬ 


tribution reported in (13). This can be seen using da m = fmnr^sinddddLp = —f nm r^ m dVL r 


and arbitrary choosing p n = zp n . With these choices, we thus obtain 
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FIG. 2. Geometry used to evaluate the integral of Eq. ©• On the left, the differential vectorial 
area da m is shown which as usual points to the outside of the closed surface Sq, centered at the 
position of the n-th atom. On the right, the polar coordinate system is shown, with the dipole 
moment p n chosen to be in the ^-direction. Note that the unit vector pointing to the outside is 

r mn = Aim- 


The second term of (Al) is easily seen to give a zero contribution to the integral: 

[ d 3 r m E rad , m (r n ) = -p n x [ d 3 r m (\7 m x = -p n x f da m x = 0 

\ ^ nm J J So—cil 


’ So-til 


'So—cil 


nm 
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Finally, the third term of (Al) is not contributing since the origin is not included in the 
integration domain: 


I d 1"rndd'rad,m{j'n ) Pn I 
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n i d 3 r m 5(r nm ) = 0 (A4) 

nm / J So—cil 


This demonstrates the result reported in Eq. (13). 
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